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ABSTRACT 



We report on the gamma-ray observations of giant molecular clouds Orion A and B 
with the Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Telescope. 
The gamma-ray emission in the energy band between ~ 100 MeV and ~ 100 GeV is 
predicted to trace the gas mass distribution in the clouds through nuclear interactions 
between the Galactic cosmic rays (CRs) and interstellar gas. The gamma-ray production 
cross-section for the nuclear interaction is known to ~ 10% precision which makes the 
LAT a powerful tool to measure the gas mass column density distribution of molecular 
clouds for a known CR intensity. We present here such distributions for Orion A and B, 
and correlate them with those of the velocity integrated CO intensity (Wco) at a 1° x 1° 
pixel level. The correlation is found to be linear over a Wco range of ~ 10 fold when 
divided in 3 regions , suggesting penetration of nuclear CRs to most of the cloud volumes. 
The Wccrto-mass conversion factor, Xco, is found to be ~ 2.3x 10 20 cm _2 (K km s -1 ) -1 
for the high-longitude part of Orion A (/ > 212°), ~ 1.7 times higher than ~ 1.3 x 10 20 
found for the rest of Orion A and B. We interpret the apparent high Xqo in the high- 
longitude region of Orion A in the light of recent works proposing a non-linear relation 
between H2 and CO densities in the diffuse molecular gas. Wco decreases faster than 
the H2 column density in the region making the gas "darker" to Wco 

Subject headings: molecular clouds: general — molecular clouds: individual(Orion A, 
Orion B) 



1. Introduction 



The Orion A and B clouds are the archetypes of local giant molecular clouds (GMC s) where 
interstellar gas condenses and stars are formed (e.g., iBergin &: Tafallal 120071 : iBallyl 120081 . and ref- 
erences therein). The clouds have been studied in various wavebands including millimeter obser- 
vations of the transition lines between CO rotational states, especially from J = 1 to J = (e.g., 
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2011 ). infrared emission (e.g.. iBeichman et al.lll988l). attenuat i on of star light (e.g.. 



Sanders et alJll984 : Maddalena et alJll986l:lDame et alJll98l l200ll : IWilson et al.ll2005l: iFukui et al 



Dobashi et al 



2010; 



Dobashi 



2005), and near infrared extinction ([Rowles Froebrichl |2009| ; iFroebrich &: Rowles 
201 lh . The two clouds are prime targets for the Large Area Telescope (LAT) on-board the Fermi 



Gamma-ray Space Telescope (Fermi) in the research of molecular clouds and CR interaction be- 
cause they lie is olated from the Galactic plane and no intense gamma-ray point source overlaps 
with the clouds (|Abdo et al.ll2009dl : lAbdo et al. IhoiObl ). 



Gamma rays from the Orion-M onoceros region wer e first detected by C OS-B in the energy 
range between 100 MeV and 5 GeV (|Caraveo et al.lll98ol: iB loemen et al.lll984l). EGRET detected 
gamma rays in the range between 100 MeV and ~ 10 GeV (jDigel et al.l Il995l . Il999l ). In these 
studies, the gamma-ray intensity distribution in a region including Orion A, B and Monoceros R2 
was fitted with three independent contributions, one proportional to the atomic hydrogen (H I) 
column density, another proportional to the CO line intensity (Wcofl and the last, a presumed 
isotropic distribution. Under the assumptions that Wqo traces the H2 column density, the CR 
spectrum doesn't change in the region and H I spin temperature (Ts) is constant, the ratio Xqo 
was determinecH, from the ratio of the gamma-ray intensities associated with the H I and CO 



distributions, to be X C o = (2.6 ± 1.2) x 10 20 cm~ 2 (K km s" 1 )" 1 (|Bloemen et al.lll984T ) and X co 



(1.35 ± 0.15) x 10 20 cm 2 (K km s 1 ) 1 (jDigel et al.lll999l ). The ratio was not separately measured 



for the three clouds, Orion A, B and Mon oceros R2, due to t he limited statistics and spatial 
resolution of the instruments. We note that IStrong et al.l (jl988l ) determined Xqq ° n the diffuse 



Galactic gamma r ays observed by COS-B to be X CO = (2.3 ± 0.3) x 10 20 cm" 2 (K km s" 1 )" 1 and 



Dame et al.1 (|200ll ). by comparing smoothed infrared intensity and Wqq distributions across the 



Galaxy, determined it to be X C o = (1-8 ± 0.3) x 10 20 cm _2 (K km s 



-In— 1 



Since the publications on the EGRET data (jDigel et al 



1995 



1999j), much progress has been 



made in studies on Orion A and B: new observational data became available (e.g..lDame et al.ll2 00 1 ; 
Lombardi k Alvesll200ll:lwilson et all2005l:lKalberla et al.ll2005l : lDobashi et al.ll2005l : lRowles & Froebrich 



200 



9: 



Froebrich &: 



Rowles 



20101 : lDobashill201ll ): study of the molecular clouds was renewed (e.g., 



Wilson et alJ^OOsT Ballyl 20081 ) : anew modeling of the Galactic diffuse gamma- rav emission was pro- 
posed incorporating large-scale CR propagation (IStrong Moskalenkdll998l : IStrong et al.l 120 00); 



theoretical calculations of collisional CO rotational-level excitation were revisited (IMengel et al 



1 We define Wco as the velocity-integrated intensity of the transition line between J = 1 to J = in 12 C le O. 

2 Our Xco is a factor converting Wco to mass column density measured in units of the proton mass in cloud 
concentrations predominantly consisting of H2. In some literature Xco is used as the factor converting V^co to 
H2 column density. Where Wco traces H2 accurately and the chemical state of hydrogen is predominantly in H2, 
the 2 definitions are expected to agree. The helium and heavier atoms are assumed to be mixed uniformly in the 
interstellar gas with the solar abundance. We warn readers that comparison of Xco values calculated on different 
CO surv eys and gamma-ray o bservations are not straightforward due to differences in their calibration procedure 



(e.g. see Bronfman et al. 19881 . for the CO calibration factor) as well as in the assumptions on the CR composition 



and the associated cross-sections. 



2001 



2007; see also Kalberla et al 



Flower 2001; Cecchi-Pestellini et al 



2005; 



Liszt 



2002 



Balakrishnan et al.ll20ol ; Iwernli et aOhood ; Ishepler et al 
200G, 2007) and the distance to the Orion nebula in the 



Orion A cloud was m easured accurately (jSandstrom et al.l 120071 ; iMenten et al.l 120071 ; iHirota et al 
2007l : lKim et al.1120081 ). 



The Fermi Gamma-ray Space Telescope mission, launched on 2008 June 11, has been sur- 
veying the sky with the Large Area Telescope (LAT) since 2008 August. Its wide field of view, 
large effective area, improved spatial resolutio n, and broad energy coverage provide much higher 
sensitivity relative to its predecessor EGRET (jAtwood et al.1 12009 : lAbdo et alj|2009al ). 



Studies based on EGRET observations have established that gamma rays from Galactic molec- 
ular clouds are dominated by neutral pion decays (which we refer t o as the "pionic gamma rays" or 



1993 



Digel et al 



"pionic emi ssion" ) in the energy band between 0.2 GeV and 10 GeV (jBertsch et al.l 
1995 . 1999 ). Orion A and B are located far (~ 8.8 kpc) from the Galactic centej§ and displaced 
from the Galactic plane by ~ 140 pc. The two clouds are only ~ 400 pc away from the solar 
system where spectra of CR species upto the sub-TeV domain are predicted to be similar to those 
measured directly at the Earth after correction for the solar modulation. 

We can now analyze Orion A and B through the high-energy gamma rays detected by the Fermi 
LAT in the light of the recent developments and study the relation between Wco an d mass column 
density (or Xqo) 111 various parts of the Galaxy and obtain the total mass of the cloudfl The 
improved spatial resolution and higher gamma-ray statistics provided by the Fermi-LAT allow us 
to determine the relation on angular scales of 1 x 1 deg 2 (pixels), without being directly affected by 
the thermodynamical, chemical, or radiation environment inside the Orion clouds, albeit within the 
limited angular resolution of the Fermi LAT and uncertainties due to any unresolved weak sources 
and CR flux variation. The results can be used conversely to study various environmental effects 
on Xqo in the translucent parts of clouds where most gas in Orion A and B resides and where the 
x r,a 



1988 



factor has not been s t raightforward to deri v e fe.g.. Ivan Dishoeck & Black 



Bolatto et aD 1 19991 : Mag nani et al 



20071 : iBurgh et all 120071 : IWall 



200 



7: 



2003; 



Bell et al 



Sheffer et all 



20081 ) . 



1986; 



Magnani et al 



20061 : ISnow fc McCalll boOfil : bell et al 



Theoretical analyses have long su ggested that Xgo depends on the environment and the Wco 



N{B. 



1988 



Taylor et al 



1993 



2 ) relation may be nonlinear (e.g., Kutner &; L eung 1985; Dick man et al.lll986 ; Malonev Sz B lack 



Bolatto et al 



1999 ; Magnani et al 



2003 



Bell et al 



200 



7: 



Burgh et al 



20071 ). Sugge stions have also been made that Xgo depends on the relative abunda nces of CO, C I, 



and C II (e.g. Ivan Dishoeck fc Black 



1988 



Hollenbach et al. 



1991 



Kopp et al.ll2000l ) . The existence 



of gas not traced by H I and CO at the interface between the two phases (the "dark gas") has been 



3 We assume the distance between the Sun and the Galactic center to be 8.5 kpc and the Galactic rotation velocity 
near the Sun to be 220 km s . 

4 The mass of Orion A and B is distributed mostly in the column density range corresponding to a "translucent" 
cloud whose line-of-sight visual attenuation (Ay) is t ypically between 1 and 5 mag and has n(H2) typically between 



100 and 2000 cm -3 (e.g. Ivan Dishoeck fc Black| [l988i. 



discovered ([Grenier et al.ll2005l ; lAde et al.ll201ll ). The relation between the fraction of carbon in CO 
and H2 density in translucent and diffuse clouds has b e en updated ba s ed on observations and nu - 
merical simulations, for example, bv lBurgh et al.l (|20ld ); IWolfire et al.l (|2010l ); lGlover et al.l (|20ld ). 
Our results will be interpreted in the light of these recent works. The Wqo-N(H.2) relation will be 
characterized including the "dark gas," and the measured mass column density will be related to 
the Ay value at which the relation is predicted to become non-linear. 



In this paper we analyze diffuse gamma rays spatially associated with the molecular cloudqfl 
Orion A and B, extract their pionic gamma-ray com ponents, obtain ma ss dis tributions, and com pare 
them with those predicted for Wqo measured by iFukui et al.l (|201ll ) and iDame et al.l (|2001l ) . In 
Section[2]we describe the gamma-ray event selection applied in this analysis. The analysis procedure 
is described in Section[3]in 4 subsections: the spatial templates used to extract mass column density 
associated with multiple emission components are given in Subsection 13. lj energy-binned spatial 
fits on the templates are described in Subsection 13.21 the pionic emission is extracted from the 
spectra obtained in the spatial fits and Xqq is calculated thereon in Subsection 13.31 and the total 
H2 masses of Orion A and B are estimated in Subsection 13.41 In Section EJ we assess systematic 
uncertai n ties in the analyses; check the Xqq results with recent infrared excess emission maps by 
Dobashil (|201ll ); summarize the results; and interpret them in the light of recent studies of the 



relation between the H2 and CO fraction in the translucent clouds. The paper is concluded in 
Section [5j 



2. Observations and Data 

The data used in this analysis were obtained in the nominal all-sky survey mode between 
2008 August 4 and 2010 March ifl W e select events classified as Pass6 D iff use class which has a 



high gamma-ray purity (jAtwood et al.ll2009l ). Among the events, we limit the reconstructed zenith 
angle to be less than 105° to greatly reduce gamma rays coming from the limb of the Earth's 
atmosphere. We select the good time intervals (GTIs) of the observations by excluding events that 
were taken while the instrument rocking angle was larger than 52°. Another cut is made on the 
reconstructed gamma-ray energy at E m i n = 178 MeV and E max = 100 GeV to reduce systematic 
uncertainty of the LAT effective area and residual background events induced by CRs. Gamma 
rays in a rectangular region of 30° x 30° centered at (£ = 210°, b = —20°) are then selected for 
later analyses. We refer to the region as the region-of-interest (ROI) and the set of events as the 
data set. 



5 By molecular clouds we mean spatially identified clouds without distinguishing the small admixture of atomic 
and ionized hydrogens therein. 

6 Mission Elapsed Time 239,557,413 s through 290,000,000 s where zero is set at 00:00 UTC on 2001 January 1. 
During the period, the LAT was operated in the survey mode with the rocking angle 35 deg (2008 August 4 to 2009 
July 9), 39 deg (2009 July 9 to 2009 September 3) and 50 deg (2009 September 3 to 2010 March 11). 



The data set consists of 1,132,436 events of which 901,929 are between 178 MeV and 1 GeV, 
224,753 between 1 GeV and 10 GeV, and 5,754 between 10 GeV and 100 GeV. They are binned in 
150 x 150 equal-area pixels (Hammer-Aitoff projection) in Galactic coordinates with 0.2° gridding 
on their reconstructed arrival directions, and in 22 logarithmic bins between E m [ Q = 178 MeV and 
-Etaax = 100 GeV on their reconstructed energies. 

The map of counts integrated over the energy range of the data set is shown in Fig. [TJ We can 
visually identify Orion A and B near the center of the region and the outer Galactic plane in the 
upper part. We note that Monoceros R2 is also visible between Orion A/B and the outer Galactic 
plane. 

3. Analyses 

The analyses presented here begin by finding the relationship between the spatial distributions 
of gamma rays and Wqo , the most widely used proxy of H2 , in the Orion clouds and by studying the 
proportionality between the two and its spatial dependence within the Orion clouds. The analyses 
proceed in 3 steps. 

In the first step, the spatial distribution of the "background" gamma rays, i.e., the gamma 
rays not associated with the H2 clouds, is determined by using spatial distribution templates, for 
the H I gas, for the inverse Compton scattering (IC) component, for the point sources, and for a 
presumed isotropic component (Subsections 13. ip . We then fit, in Subsection 13.21 the gamma-ray 
spatial distribution in each of the 22 energy bins as a sum of the "background" distribution and a 
distribution tentatively associated with the H2 gas (H2-template). The "background" is subtracted 
from the measured gamma-ray intensity distribution and the remainder is defined as the gamma-ray 
intensity distribution associated with the H2 gas with which Wco is correlated pixel-by-pixel. We 
note that the gamma-ray intensity measures the mass column density in the H2 gas for a known 
CR spectrum. We repeat the fit with 2 alternative H2-templates. 

In the second step (Subsections 13.31 and 13.41) . the energy-binned gamma-ray emissivity for the 
H2 gas (Bi in eq. ([I])) are assembled as the gamma-ray spectrum for each of the 3 H2-templates. 
The spectrum is then fitted as a sum of the gamma rays produced in the pionic and bremsstrahlung 
processes. 

In the third step, the gamma-ray intensity distribution associated with the pionic emission is 
converted to the mass column density. The Wco- m ass conversion factor (Xco) is calculated via 
two methods, one by comparing the gamma-ray counts associated with the H I gas and with the H2 
gas (the H2/H I method) and the other by dividing the gamma-ray counts of the pionic emission 
by the number of pionic gamma rays expected per unit gas mass (the pionic method). In the 
first method, we assume the CR spectrum is uniform in the local H I region within Galactocentric 
radius of 8 — 10 kpc (see Subsection I3.1.1j) and in the Orion clouds. In the latter method, we 
assume the CR spectrum including its absolute flux is known in the Orion clouds. We validate 



these assumption using GALPROP. 



We use GALPROP (jStrong fc Moskalenko 



1998 



Strong et al.l |2000| ) with the parameter set 



labeled as GALDEF 54_77Xvarh7S. This parameter set is kn own to reproduce r easonably well the 
Galactic diffuse gamma-ray emission observed with the LAT (kbdo et al.ll2()()9lfl. We refer to the 
results obtained by running GALPROP with this parameter set as the GALPROP results in this 
paper. 



3.1. Spatial distribution templates 

Initially we assume the gamma-ray emission from the ROI to be made of 4 "background" 
components and one "signal" H2 component, each emitting gamma rays with a characteristic spatial 
distribution. The 4 "background" components are spatially associated with the diffuse H I gas, the 
inverse Compton (IC) scattering by electrons^! off interstellar radiation fields, the point sources, and 
the sum of extragalactic diffuse emission (including unresolved sources) and backgrounds induced 
by CRs in the instrument. We assume the last sum to be isotropic and refer to it as the isotropic 
component. We ignore the contribution from ionized hydrogen gas (H II) because its dens ity is low 



when averaged in 1 x 1 deg 2 pixels (< 0.5cm 3 ) and its total mass is negligible in the ROI ( Gordon 



1969 



O'Dell 



2001 



All spatial components except for the IC component are assumed to have, individually, an 
energy-independent underlying spatial distribution in Galactic coordinates (l,b). Another impor- 
tant underlying assumption is that the nuclear CR spectrum is uniform over the ROI. We make 
spatial templates for the 22 energy bins by convolving the spatial distributions with the energy- 
dependent point spread function (PSF) and exposure for the individual energy bins. Hence the 
spatial templates are energy dependent. In Subsection 13.31 we will show that the spectra of the 
gamma-ray emissions associated with the H I and H2 gas consist of the pionic and bremsstrahlung 
components. 

The gamma-ray intensity I~(l,b) for the i-ih energy bin is interpreted as the sum of the five 
contributions, each being the product of the normalization factor for the i-th energy bin and the 
spatial template. 

J 7>i (i, b) = AiN(R b) + B, t N(R 2 )(l, b) + ICS, b) + £(Q;<W + A, (1) 

3 

The normalization factors are: A4 for the H I gas; Bi for the component associated with clouds 
consisting predominantly of H2; ICj for the inverse Compton component; Cij5i j: bj for the j-th point 



7 A detailed description of GALDEF files can be found at http://galprop.stanford.edu 
8 We refer to electrons as a sum of e + and e~ . 



source at bj); and D{ for the isotropic component which is assumed not to depend on (l,b). The 
normalization factors are determined independently for the 22 energy bins. We note that ICj are 
fixed at the values given by GALPROP, because the spatial distribution is highly correlated with 
the isotropic component, and the IC component is sub-dominant in the ROI. 



Later in Subsection 13.21 we will explore 3 templates for H2 , two based on Wqo and one on 
Wqo phis the "dark gas" proposed bv lGrenier et al.1 ( 20051 ). 



3.1.1. Diffuse H I gas template 



Atomic hydrogen gas (H I) is bro adly distributed in the Galaxy with a total mass exceeding 
that of molecular hydrogen (H2) (e.g.. lFerrierdl200ll ; ISnow &: McCalll 120061 ) . In the outer Galaxy 
where the Ori on clouds are located , the m ass column density of H I is lower than that of H2 at the 
Orion clouds (jKalberla et alJboid . 120051 ) 



We used the Leiden/ Argentine/Bonn (LAB) survey data (jKalberla et al.ll2005l ) corrected for 
optical thickness by adopting a constant spin temperature (Tg) of 125 K as the H I gas spatial 
distribution template (see Fig. [2^,) . The LAB intensity distribution is divided into five annuli cen 



tered at the Galactic center as has been done in other Fermi diffuse emission analyses (jAbdo et al 



2010a ). Their inner and outer Galactocentric radii (R) are: 8 to 10, 10 to 11.5, 11.5 to 16.5, 16.5 
to 19, and 19 to 50 kpc. The line-of-sight velocity distribution of the H I gas in the Orion region 
overlaps that of the CO gas associated with the Orion clouds and that of the local H I annulus 
(R = 8 - 10 kpc) quite well. 

Gamma-ray contributions from all the H I annuli overlapping our ROI have been included 
in the analyses. In the fitting, the CR intensity is treated independently at each annulus. The 
contributions from annuli other than the local one (R = 8 — 10 kpc) are through the periphery 
of the LAT PSF and less than ~ 5% in gamma-ray counts. Hence our analyses are insensitive to 
variation in the CR intensity and/or spectrum among the neighboring annuli. 

The spin temperature of H I gas, Tg, is not well constrained in the region nor known to be 
unifo rm over the ROL it s quoted value in the literature ranges between ~ 90 K and ~ 400 K 
(e.g. iMohan et all l2004al Jbh. We estimate, later in this paper, the contribution to the overall 
systematic error from this uncertainty by repeating the analysis for Ts = 250 K and 90 K. No 
significant concentrat ion of cold H I is known around Orion A and B at large-scale (> 1 x 1 deg 2 ) 
(jKalberla et al.ll2010l ). An exploratory study o f cold H I mixed i n selected H2 cloud cores has found 
the mean HI fraction to be less than 0.5% (jKrco et al.1 120081 ). So we can safely ignore such a 
mixture in the analysis. 



Gamma rays are produced in the H I gas through the pionic and bremsstrahlung processes 
with intensities proportional to the CR nuclear and electron spectra in the gas, respectively. 



3.1.2. Molecular cloud template 



We try 3 H2 templates to represent the H2 spatial distribution in the ROI. In making the 
templates we assume that the H2 col umn density is pro portional to Wco measured by two CO (J = 
1 — >■ 0) surveys, one from NANTEN (jFukui et al J 1 20 111 ) covering the ar eas around the Ori on clouds 



with effective resolution of 4' and the other being the Galactic survey by lDame et al.l (|200ll ) covering 
the ROI with angular resolution of 8.7'. The spatial distributions indicated by the two surveys are 
mutually consistent at the angular scale of the LAT PSF except for the overall normalization. 

The first H2 template, H2-template-l, is made by combining the two surveys and accounting 
for their relative intensity sca l es (Fi g. [2)3) : NANTEN Wqo for the area defined by the solid white 
line and that by I Dame et al.l (|200ll ) for the rest of the region. We refer to the 3 regions defined 
by dashed lines in Fig. [3)3 as "the 3 Orion regions" hereafter^. We scale the NANT EN data by 
the factor 1/1.11 to adjust t he intensities to a common scale with lDame et al.l (|200ll ) because the 
survey by lDame et al.1 (|200ll ) has been widely used in gamma-ray analyses. 



We first start the analyses by setting one common B{ factor for Wqq in the ROI (H2-template- 
1), or equivalently, one common A"co for the entire ROI. In the second H2 template, H2-template-2, 
the Wco distribution is divided into 4 regions (the 3 Orion regions and the rest of the ROI) and 
allow Bj, or equiva lently Xqo, to be different in each region. We add a "dark gas" template 
(jGrenier et al.ll2005l ) to H2-template-l to make the third H2 template, H2-template-3 (Fig. [2^ and 
Et). The normalization is set free for the 2 templates. 



These spatial templates are described further in the subsections to follow. 



3.1.3. Inverse Compton template 



The inverse Compton component is known to be minor around the Orion clouds. We use the IC 
spect rum and spatial dis tribution given by GALPROP where the interstellar photon fields are taken 
from lPorter et aL ( 2008 ). The typical Galactic-scale IC intensity in the region is ~ 5 times smaller 
than the isotropic component described later, and their spatial and spectral distributions are similar 
in this region. Possible local enhance ment is the IC emi ssion around the Orion Nebula (M42) 
where strong ultravio let emission (e.g. iMurthy et al.l 120051 ) and moderate infrared emission (e.g. 
Prisinzano et al.ll2008l ) exists. Ac cording to our calculatio n, such IC emissions are not detectable 
with the current LAT sensitivity (jOrlando &: Strong| |2008). 



9 The boundaries are: Orion A Region I (217° > t > 212°, -23° <b< -16°), Orion A Region II (212° > I > 205°, 
-23° < b < -16°, excluding the overlap with Orion B), and Orion B (209° >£> 203°, -18° < b < -13°) 



3.1.4- Point sources in the Orion region 



More than 1400 point sources are reported in the First Fermi LAT Catalog (lAbdo et al. 



2010H ). Among them, 30 point sources are in our ROI, (l,b) = (210 ± 15°, -20 ± 15°). There 



are an additional 29 sources within 5 deg of the ROI. In the likelihood fit to be discussed later, 
the normalization is set free, energy-bin by energy-bin, for 25 high-confidence sources in the 
ROI; the indexes and nor malizations are fixed to the values given in the First Fermi LAT Cat- 
alog (jAbdo et al. I l2010bl ) for those outside of the region. There are 5 low-confidence sources 
(or candidates) overlapping with the clouds: they ar«0 1FGL J0540.4-0737c, J0536.2-0607c, 
J0534.7-0531c, J0541.9-0204c, and J0547.0+0020c. Their fluxes are all low and labeled as "c" 
in the catalog, meaning either their flux estimates are uncertain, or they can be artifacts resulting 
from incorrect modeling of the Galactic diffuse emission. We fit the spatial templates and analyze 
the spectra in the 3 Orion regions with and without them. The results we quote will be for the 
analyses without them: we include their possible contribution in the systematic error. 



3.1.5. Isotropic component 



In the present analyses, the extragalactic emission a nd residual CR background in the data are 
not s e parated but treated as a single isotropic component ( Abdo et all2009c . 2010c : Ackermann et al 



2010 



Abdo et al.ll2010al ). The total flux of the component at 1 GeV is ~ 25% of that associated 
with H2 when averaged over the 3 Orion regions (subtending ~ 30 msr) defined around Orion A 
and B (see Fig. [3b). 

The residual background in the Pass6 Diffuse class consists of CR-induced events misclassified 
as gamma rays and CRs that converted in the passive material just outside of the LAT without 
leaving a signal in the anti-coincidence detector (lAtwood et al.l 12009). When averaged over many 
orbits of observations, the residual background can be approximated as isotropic. 



3.2. Fit to the Spatial Distribution 

All spatial templates described in the previous subsection were convolved with the LAT ex- 
posure and PSF. The spatial fit is made using the binned likelihood program gtlike included in 
the Fermi ScienceToola 11 ! and the 4 normalizations (Ai, Bi, Cy and D~) in Eq. (pQ) are determined 
independently for the 22 energy bins. We note again that ICi are fixed at the values given by 
GALPROP. Each H I annulus has a separate Aj. We report only Ai for the local annulus as others 
are not determined well because they lie mostly outside of our ROI. 



10 No new sources have been added in this region in the Second Fermi LAT Catalog 1 Abdo et al, 
11 We use ScienceTools version v9r!6p0 with P6-V3-DIFFUSE instrument response functions. 



2011 



Our scientific interest is to study the contributions from the gas concentrations identified as 
Orion A and B, which are believed to be predominantly H2. We consider, hence, the sum of the 
HI, IC, point-source, and isotropic components as the "background" which is determined by fitting 
the observed gamma-ray distribution for each of the 22 energy bins. In the fits, we assume that H2- 
template-1, or the Wco distribution, represents approximately the H2 distribution. The gamma-ray 
distribution associated with the H2 gas can be extracted less dependently on yet-unknown H^-Wco 
relation by subtracting the "background" from the observed gamma-ray distribution. 

We define 2 improved H2-templates, H2-template-2 and 3 after the initial analysis on H2- 
template-1. The spatial distribution is not proportional to Wco for the 2 improved templates and 
hence the "background" is different for each H2-template by a small amount. The difference is 
however negligible. 

3.2.1. Spatial Fit with Wco of One Xco: ^-template- 1 

We use H2-template-l as an approximation for the H2 gas distribution and fit Eq. (pQ) to 
determine the "background". The energy-summed gamma-ray distribution after subtracting the 
"background" is shown in Fig. [3^i and that of the Wco-based model, or the product of Y>Bi in 
Eq. dl]) and H2-template-l, is given in Fig. [3)3. The two count distributions are correlated pixel- by- 
pixel (lxl deg 2 ) in the 3 Orion regions in Fig. 0^. We expect a good linear correlation between 
the two if Wqo is a good tracer of H2. 

We note first that the correlation is fairly linear and gives a correlation coefficient!* 2 ! of 0.93. 
We then note that the correlation significantly improves if we separate the Orion clouds into the 
3 Orion regions, Orion A Region I (black solid line) and II (red dashed line), and Orion B (blue 
dotted line). The correlation coefficients for the 3 Orion regions are 0.98, 0.96, and 0.98, and the 
best-fit slopes are 0.72, 0.99, and 1.25, respectively. 

The large difference (~ 40 — 60%) in the best-fit slope suggests that the mass column density 
in Orion A and B cannot be simply derived using the same value of -Xco- We find more gamma 
rays in Orion A Region I per Wqo than in Orion A Region II and Orion B, suggesting Xco is 
different in the 3 Orion regions, or that some fraction of the H2 gas is not traced by Wqo provided 
a uniform CR density. We explore these two possibilities by redefining the H2-template. 

3.2.2. Spatial Fit with Wqo °f 4 different Xqo values: ^-template- 2 

Based on the relation found between the spatial distributions of the gamma-ray intensity 
associated with the H2 gas and the Wco-based model (H2-template-l), we make a second template, 



The correlation coefficient is defined as T,(x — x)(y — y)/^Y,(x — :r) 2 E(y — y) 2 



H2-template-2, that will delineate the H2 column density more faithfully. In the template we divide 
the ROI into 4 regions, the 3 Orion regions and the rest of the ROI, and allow Bi to be different in 
each region, or introduce 4 B^s. 

The fitted results for A{ (H I) and Bi (H2-template-2) in Eq. ([TJ are listed in Table Q] after 
combining the highest 10 energy bins into 3 bins. The gamma-ray count map is shown in Fig. [3h 
is the sum of the 4 B^s multiplied with the corresponding components of H2-template-2. We note 
that the 3 Orion regions mix to some degree through the Fermi PSF. The correlation between the 
gamma-ray distribution associated with H2 and the H2 template improved as shown in Fig. Hb: the 
best-fit slopes for Orion A Region I, Region II and Orion B are 0.95, 0.94, and 1.03, respectively, 
while the correlation factors remain almost the same, 0.98, 0.99, and 0.96, respectively. 

The Xqo for the 4 regions can be calculated directly as the ratio of Bi to 2Ai (the H2/H I 
method) or by extracting the gamma-ray emission in the regions (the pionic method). The results 
from the former are given in Table [2] together with those from the latter which will be described in 
Subsection 13.31 



3.2.3. Spatial fit with Wco and "dark gas": H2-template-3 



Grenier et al.1 (J2005J) found that a significant fraction of local diffuse gamma-ray emission 
observed by EGRET is not associated with either H I or VFco, but rather with the dust map 
traced by thermal infrared emission given by lSchlegel et al.1 (119981 ). The missing gas component is 
often referred to as th e "dark gas". Other LAT observations ha ve found gamma rays associated 
with such "dark gas" (jAbdo et al.1 l2010a : lAckermann et al.l l2010l ) . We note recent measurements 
of attenuation or reddening of background stars have also detected gas concentrations not traced 



well by Wqo (jDobashi et al. 



2005 



Rowles Froebrich 



2009 



Dobashi 



2011 



Ade et al 



2011 



)• 



We make a third template, H2-template-3, that can bring out the true gas distribution associ- 
ated with the Orion clouds and enhance our understanding of the VFco-to-H2 relation by introducing 
the "dark gas". The new H2 template consists of H2-template-l, or Wco, and a "dark gas" spatial 
template with a normalization factor for each. 



Our "dark gas" template has been produced following the prescription given by I Grenier et al 

(|2005l ) and referred to as E(B — V) res . It is a residual map obtained by subtracting the b est- 
fit linear co mbination of iV(H I) an d Wco from the E{B — V) map of Schlegel et al. (jl998l ) as 
described in lAckermann et al.l (j2010l ). Fig. [2fc shows the E(B — V) i e s map around our RO I. There 
is a problem with the color temperature correction of the map by I Schlegel et al.l (| 19981 ) around 
the OB associations in the Orion A and B clouds, and thus E(B — V) Tes value is negative in these 
points. We masked out these pixels in the E(B — V) Tes map by setting the corresponding values 
to zero. 



The results for Ai (H I) and Bi (2 normalizations, one for Wco an d one for the "dark gas" ) in 



Eq. dU) are listed after combining the highest 10 energy bins into 3 bins in Table O The distribution 
of the gamma-ray counts associated with H2-template-3, the sum of the counts associated with Wco 
and the "dark gas", is given in Fig. [3jd. The correlation between the extracted gamma-ray counts 
and the model counts improves as shown in Fig. 0b, bringing the correlation coefficients to 0.99, 0.99, 
0.97, and 0.98, for Orion A Region I, Region II, Orion B and the sum of the 3 regions, respectively. 
The improvement in the correlation, or equivalently in the spatial fit, comes from inclusion of 
E(B — V) ies which has the largest contribution in the Orion A Region I seen in Fig. [2b. 

3.2.4- Summary of the Spatial Fits 

The relative likelihoods of the spatial fits with Eq. (fTJ) in the ROI are compared among the 3 
H2-templates in Fig. [5] for the 22 energy bins. The "dark gas" template (H2-template-3) gives the 
best fit in almost all energy bins and the 3-Xqo template (H2-template-2) gives the second best 
result. The improvements relative to H2-template-l are statistically significant. 

The residuals of the fits with the 3 templates in the ROI are given in Fig. [6l The rectangular 
boundaries of the 3 Orion regions shown in Figs. [3b, c, d are replicated in the figure. The residuals 
are significant within the Orion regions for H2-template-l (Fig. [6b) but not for the other 2 templates 
(Fig. [6b, c), which is consistent with the improvement we saw in Fig. [H The difference in the 
residuals for H2-template-2 and H2-template-3 in the Orion regions is not significant relative to 
the systematic uncertainty discussed in the next subsection. We find that the large improvement 
H2-template-3 has brought relative to H2-template-2 in Fig. [5] comes primarily from outside of the 
3 Orion regions, especially in the Monoceros R2 region and in the northern region adjacent to the 
Orion B: the template adds "dark gas" in that part whereas the other templates only modify the 
3 Orion regions. 

The value of Xqo has been calculated by the H2/H I method by taking the ratio of Bi to 
2Ai for the parts associated with Wco i n the ^-templates and listed in Table (2[ In the pionic 
method of evaluating Xqoi however, the pionic component must be extracted out of the gamma-ray 
spectrum associated with the H2-template as will be described in Subsection 13.31 We will discuss 
the systematic errors in evaluating Xqo an d possible interpretations of the results in Section [U 

3.3. Analyses of Spectra 

The spectra associated with the H I and H2-template-l, with the H I and H2-template-2, and 
with the H I and H2-template-3 are obtained by assembling the fitted results for the respective 
templates, Aj and Bi, as shown in Figs. [HJ and[9l respectively. The spectra are fitted as a sum 
of the pionic and bremsstrahlung components. The gamma-ray spectra associated with the spatial 
templates (H I, inverse Compton, isotropic, and sum of Xco x Wco) are plotted for the 3 Orion 
regions in Fig. [TUb and b. We analyze for the gamma rays associated with the 3 H2-templates in 



this subsection. 



3.3.1. Fit with Gamma-ray Emission Models 



The spectral template of pionic gamma rays has been calcu 
ray inclusive cross-section for p — p interaction parameterized by 



ated by convolving the gamma- 



Kamae et al.1 (120061 ) and the CR 



The proton flux is predicted in 



proton spectrum predicted by GALPROP at the Orion clouds 
the Orion clouds (R = 8.8 kpc, Z = —0.14 kpc) to be ~ 8% smaller than that at the solar system 
(R = 8.5 kpc, Z = 0.0 kpc) where the GALPROP proton spectrum has been determined by the CR 
data taken at the Earth. The value at the Orion cloud s is consistent with that determined using 
the gamma rays from the local H I (jAbdo et al.ll2009d ). The good fit to the data seen Figs. [H 
and [9] supports GALPROP's prediction of CR spectral shape in the Orion region and the overall 
modeling of Eq. ([TJ . 

Bremsstrahlung emission induced by CR elec trons interacting with gas is calculated in GAL 



PROP using recent bremsstrahlung calculations ([Strong &: Moskalenko 



1998 



Strong et al.l 1200(1 . 



and references therein). The electron injection spectrum in our GALPROP calculation had been 
adjusted to repr oduce, approximat ely, the power-law index of the electron spectrum measured by 
the Fermi LAT (JAbdo et al.ll2009d ). In addition, the normalization of the spectrum is adjusted to 
reproduce the LAT observed gamma-ray flux at a low-energy band. In the spectral fits described 
below, we kept the electron-to-proton ratio, or equivalently the bremsstrahlung-to-pion ratio, fixed 
to the value given in GALPROP. When we refer to the gamma-ray emissivity per atom or molecule, 
we do not differentiate the underlying processes, but rather the sum of the bremsstrahlung and 
pionic contributions. 

The spectral fit of the H I component is reasonable for all 3 H2 templates (x 2 = 17.7, 9.9, and 
17.1 for /dof = 14 respectively) as shown in Figs. [Th., [Hh. and [9h. Our pionic flux ass ociated with 



H I is consistent with that obtained in the Fermi study on the local interstellar gas (jAbdo et al 
2009q ) as overlaid in Fig. [Hh.. We note however that there may be a small offset between the two 



as will be discussed later. The spectra associated with molecular clouds are also fitted well by the 
3 H2-templates as shown below. 

The mass-to- Wco ratio, Xqo, can be obtained by comparing the assumed pionic gamma-ray 
emissivity per H atom with the observed gamma-ray emissivity per Wco as shown in Figs. [TJa, 



13 In GALDEF 54_77Xvarh7S, the CR proton flux was artificially multiplied by 1.15 to reproduce gamma-ray 
observations by Fermi. The factor originates from the underestimate of gamma-ray emissivity for He and heavier 
atoms in the i nterstellar medium (ISM) in GALPROP. Instead of using the 1.15 correction factor, we combined the 
calcu l ation by iGaisser fc Schaeferl (|1992T ) for contributions from CR He and heavier atoms, and the calculation by 
Mori 1 20091 ) for heavier atoms in the ISM. Hence the total gamma-ray emissivity per H atom is 1.70 times larger 
than that for p — p collisions only. The difference between the total gamma-ray emissivity in the two literatures is 
~ 5%, which is taken into account in the systematic uncertainty. 



03-d and [9)3. The former is calculated in the unit of MeV s 1 sr 1 , and the latter is measured in 
the unit of MeV~ 1 s _1 sr" 1 (2 x 10 20 cnr 2 (K km s -1 ) -1 ) -1 . Thus, X C o/2 of the clouds is derived 
by dividing the latter by the former. 

The results of the spectral fit for the H I component are not used to determine -Xco in the 
pionic method. The fits to the spectral components shown in Figs. Ek, [Bh^ and[9h are only to check 
overall consistency of our analyses. Their normalizations are consistent within the uncertainty in 
the H I column density discussed in Section HI 

3.3.2. Spectra obtained with H2 -template- 1 

The fitted spectra are plotted as sums of pionic and bremsstrahlung emissions in Fig. [7k., b 
for the H I spatial template and the H2-template-l (Orion A Region I, II, and Orion B combined), 
giving % 2 /dof of 17.7/14 and 20.2/14, respectively. 

We give the -Xco value obtained from the fitted pionic spectra in Table [2j Since the fit is 
substantially poorer than those for H2-template-2 and 3 (see Fig. [5|), the value should be taken just 
as a reference value. For this reason we do not quote systematic errors in the table. 

3.3.3. Spectra obtained with ^-template- 2 

The fitted spectra are plotted as sums of pionic and bremsstrahlung emissions in Fig. [8)3, c, d 
for Orion A Region I, II and Orion B, giving x 2 /dof of 14.0/14, 18.5/14, and 10.6/14, respectively. 
The -Xco values obtained for the 4 regions from the fitted pionic spectra are given in Table [2j 

The coefficient -Xco is significantly higher for Orion A Region I than for other regions, con- 
sistent with the slopes obtained in Subsection 13.21 in the pixel- by-pixel correlation study. This also 
can be seen in the -Xco obtained with the H2/H I method. 

We note that the fraction of the H I component in the gamma-ray spectrum integrated in the 
3 Orion regions is comparable to that associated with Wco (see Fig. [TUb). This is because the solid 
angle subtended by the Orion molecular clouds is a small fraction of our 3 Orion regions in solid 
angle and the overall mass of atomic gas is greater. 

3.3.4. Spectra obtained with H2-template-3 

The fitted spectra integrated over the Wco and "dark gas" components are shown in Fig. [9J3, 
c. We give Xco f° r the ROI from the fitted pionic spectrum in Table [2j 

The Xco obtained in fits with the Wco can be compared with those obtained in similar 



analyses including the "dark gas" template: 2.0 x 10 (in the local arm), 1.9 x 10 (the Perseus 



arm and 0.87 x 10 20 (the Gould Belt) in the same unit as above (jAckermann et alJboid : lAbdo et al 
2010c! ) . 



The spectrum associated with the "dark gas" component is similar in shape to that associated 
with Wqo but about half as intense (Fig. fTUb). The 2 spectral energy densities (SEDs) become 
comparable in Orion A Region I as seen in Fig. [TUb . The "dark gas" dominates over Wqo 111 the 
pixels near the high-longitude end of Orion A and eventually Wqo diminishes in the pixels beyond 
them towards higher longitude. 

Our Xqo measurements given in Table [2] can be compared with those determined using 
the gamma-ray flux fro m the Orion-Mono ceros complex measured with EGRET: (1.35 ± 0.15) x 
10 20 cm _2 (K km s -1 ) -1 ( Digel et al.lll999i ). We note there were no Galactic CR propagation mod- 
els such as GALPROP nor CR measurements as precise as are available now: Xqo was determined 
by the H2/H I method and it compares well with the single Xqo value of 1.36 ± 0.02 obtained with 
the H2-template-l. 



3.4. Total masses of Orion A and B 

The distance from the Sun to the Orion nebula (M42) in side the Orion A has recently bee n 
measured by parallax to be 389 +2 i pc ( Sandstrom et al. 20071). 414 ± 7 pc ( Menten et al.1 120071 ) . 
437 ± 19 pc (jHirota et al.ll2007l ). and 419 ± 6 pc (|Kim et alJboOsi ). We adopted 400 pc as the 
distance to the Orion A and B clouds and used the total pionic gamma-ray fluxes obtained above 
to get the total masses of Orion A and B outsidJ^l. 

Mass estimation using H2-template-2: 



M A = (74.5 ± 1.3) x 10 3 M. 
M B = (33.5 ± 0.7) x 10 3 M 



400 



400 



where 



M 



400 



d 



x M 



05 



400 pc J ® 
and d is the distance to the clouds. We will discuss the systematic uncertainties in the next section. 

Mass estimation using H2-template-3: Addition of the "dark gas" changes the estimation of 
the Orion A and B masses by about 10%. 



M A ,w co = (55.1 ± 0.8) x 10 3 M 



400 



W e note that the spatial extent of Orion B defined here is significantly different from that used in I Wilson et al. 
(2005) because we are unable to separate Orion B from the complex cloud structures behind due to the broad PSF 
of the LAT. 



M B 



(27.6 ± 0.7) x 10 3 M 



100 



(36.0 ±0.5) x 10 3 M 40 o- 



The total mass of Orion A (= M A ,Wm + M A .r) ark ) is (82.7±T.l) x 10 3 M 400 . The Orio n A mass has 



2001 



been estimated by I Wilson et al.1 ( 20051 ). assuming 

X C o = 1.8xl0 20 cm" 2 (K km s -1 )" 1 (|Dame et al 



), to be Ma = 91. 7 x 10 3 M^n- The m ass has been estimated separately for Orion A Regions 1, 
2, 3, and NGC 2149 in lWilson et al.l (|2005l ). Our Orion A (Region I and II) includes their Regions 
1, 2, and 3 but overlaps only partially with NGC 2149. Considering the breadth of the PSF and the 
limited statistics of the data, we could not determine how much of NGC 2149 overlaps our Orion 
Region I. If we assume about one half of NGC 2149 is in our Orion Region I and the sys t emati c 
error introduced by this ambiguity is half of the NGC 2149 mass estimated by I Wilson et al.1 (|2005l ). 
the Orion A mass to be compared becomes Ma = (86.3 =L 5.4) x 10 3 M^qq. The Orion B region is 
more complex and such a comparison is very difficult. 



4. Discussion 

Although the Orion clouds lie away from the Galactic Plane and subtend relatively small solid 
angle, many Galactic and extragalactic sources contribute to the ROI through the large PSF of the 
Fermi-LAT . 

We have analyzed the observed data to extract the intensity associated with the molecular 
clouds, the 3 Orion regions in particular, by using the 3 H2-templates made from Wco on the 3 
different assumptions for each of the 22 energy bins. The ratio of the normalization factors for 
H I and H2 , A, /2Bi , gives the conversion factor of Wco to the mass column density, Xco (the 
H2/H I method). For this, the H I mass column density must be well understood from the radiative 
transfer of the H I line and the CR spectrum must be constant in the ROI. 

In the second method (the pionic method), Xco is determined by comparing the observed 
pionic gamma-ray intensities with those expected from the CR spectrum at the Orion clouds and the 
pionic gamma-ray production cross-section. For this, we have to know the absolute CR spectrum 
and flux, the instrument response function (IRF), and the pionic gamma-ray production cross- 
section, in particular the pionic gamma-ray contribution from metals in CR and ISM. 

In the subsections to follow, we evaluate uncertainties and possible systematic errors in the 
analyses, especially in evaluating Xco hi the 3 Orion regions. We then summarize the results 
obtained in this paper and present possible interpretations thereon. 



4.1. Possible Systematic Errors in the Analyses 



Systematic errors that affect the correlation measurements between gamma-ray intensities and 
Wqo are discussed in two categories: the first one applies commonly to the 3 Orion regions and 



the second affects the relation differently in the 3 regions. 



4-1.1. CR intensity at the Orion clouds 



Uncertainty in the fluxes and spectra of CRs, in particular those of protons, can affect in both 
categories. The Galactic CR protons that produce pions in our ener gy range remain in our Galaxy 
longer (~ 5 x 10 7 yrs) than electrons (~ 7 x 10 6 yrs) (|Lee et al.ll201ll ) and their flux variation within 
the Galaxy is believed to be predicted well by GALPROP. We note that the CR source distribution, 
the Galaxy size, and the CR diffusion coefficient are the important inputs to GALPROP. Using 
the CR spectrum measured at the Earth, we have calculated the CR spectrum in the Orion region 
for the 2 choices of the CR source distributions a nd the 3 choices of Galactic halo heights (2, 4, 
and 10 kpc) used in a GALPROP-based study by iLee et al.l (|201ll ). The CR spectrum does not 
change more than ~ 2% from the value used here as long as it is constrained to the measurements 
at the E arth and to repro duce the Galactic diffuse gamma-ray intensities measured by the Fermi 
LAT (see lLee et al.l (|201ll )). We also note that the gamma-ray spectrum from the local H I (typical 
distance < 1 kpc) is consistent with the CR proton flux being within ~ 10% of that at the Earth 
(|Abdo et al.ll2009ch . 



CRs could be accelerated in the clouds and/or prevented from penetrating into their cores by 
embedded magnetic field. We first note that the re are no strong non-thermal X-ray source nor radio 
SNR found in the clouds (jFeigelson et al.l 120021 . and references therein). Therefore no appreciable 
CR acceleration is likely to be taking place in the Orion clouds. The good linear correlation between 
Wco an d gamma-ray intensity seen in all 3 Orion regions (Fig. H^) confirms that the CRs effective 
in producing pions (kinetic energy > 1 GeV) are penetrating well inside the higher-density parts 
of the clouds. 

Based on these observations we assume that the CR flux in the Orion region is 8% lower than 
that at the Earth with possible systematic error of ±10% due mostly to disagreement among recent 
CR measurements at the Earth and solar demodulation uncertainties. 

Uncertainty in the CR flux at the Orion clouds contributes directly to the systematic error in 
the pionic method but indirectly in the H2/H I method. In the former, the absolute CR intensity 
is assumed to be known while the CR intensity is assumed to be the same in the local H I region 
and the molecular clouds in the latter. 



4-1.2. Uncertainty in the instrument response functions 



The uncertainty in the absolute calibration of the LAT effective area can also introduce error 
of the first kind. The effective areas were derived based on Monte Carlo stu dies of the LAT, checked 
against beam tests at accelerators (jAbdo et al.ll2009al : lAtwood et al.l 120091 ) . Comparisons between 



flight data and M onte Carlo studies have been made to quantify the systematic uncertainty in 
the effective area (jAbdo et al.ll2009bl ). At present, we estimate this systematic error to be 10% at 
100 MeV, 5% at 500 MeV and 20% at 20 GeV. 



The systematic error in the absolute energy scale has been estimated as +5/ — 10% (jAbdo et al 

2009eJ). We have refitted Xco after artificially shifting the energy scale by +5% and by —10%: the 
number of pionic gamma rays changes less than +1/ — 8% for all 3 Orion regions with all 3 H2 
templates. We include this possible error due to the uncertainty in the energy calibration when 
assessing the overall systematic error. 

The pionic method is affected directly by the uncertainty in the instrument response function 
while the H2/H I method is insensitive because it affects the denominator and numerator similarly. 



4-1-3. Uncertainty in the spin temperature of H I 



In converting the observed 21 cm line emission intensity (jKalberla et al,l 120051 ) to the H I 
column density, Tg was assumed to be 12 5 K. The range of Tg measured in the local H I gas varies 
broadly between 90 K and 400 K (e.g., bohan et al.l l2004al lbl. and references therein) while we 
have assumed a likely range for our ROI to be between 90 K and 250 K. 

We refitted the Fermi data in the ROI with these two extreme Tg values with H2-template 
2 and 3. We then calculated Aco by dividing Bi by 2Ai in Eq. (P), or by extracting the pion 
component in the spectra. The deviations of Xqq, from those obtained with Tg of 125 K are taken 
into account in the systematic errors given in Table [2) The large systematic errors for Xqo 011 
Bi/2Ai (Column 3) enter via 2^4, which depends on the absolute calibration of the H I gas density 
or Tg in the local H I. The pionic method uses the product of the CR intensity and pp — > 7 cross- 
section in place of 2Ai and is less directly affected by the uncertainty in H I gas density or Tg of 
the local H I, although the uncertainties can have a small indirect effect through the overall spatial 
fitting. This effect is much smaller than the overall systematic error and negligible. We note that 
there is s ome discrepancy be tween the gamma-ray spectra associated with H I in the ROI and the 
local H I (|Abdo et al.ll2009ch as seen in Fig. El 



4-1-4- Effect of overlapping point source candidates 

We have not included the 5 sources overlapping with the Orion clouds (Sec. I3.1.4P because 
they are all classified as "potentially confused with interstellar diffuse emission or perhaps spurious" 
(jAbdo et al. Il2010bl ) . To investigate their potential contribution we repeated the analysis including 
these sources with the fluxes and spectra listed in the First Fermi LAT Catalog. The fit with the 
pionic method gives the following A C o in unit of cm~ 2 (K km s -1 ) -1 : (2.29±0.05) x 10 20 for Orion 
Region I; (1.16±0.05) x 10 20 for Orion Region II; and (1.24±0.04) x 10 20 for Orion B. They are 2%, 



19%, and 8% less than those obtained without these point source candidates. In the present study, 
we assume they are artifacts and add +0/ — 2, +0/ — 19, and +0/ — 8% to the overall systematic 
error in the 3 regions. 



4-1.5. Overall error 

For the H2/H I method, the uncertainty in the H I mass density (~ 20%) due mostly to the 
uncertainty in Tg dominates the systematic error . Other contributions include the overlapping "c" 
sources (+0/ — 2, +0/ — 19, and +0/ — 8%) and variation in the CR intensity within ~ 1 kpc or 
between H I and the molecular clouds (±5%), making the total systematic errors for the 3 Orion 
regions to +25/ - 28, +25/ - 44, and +25/ - 33% as given in Column 3 of Tabled 

For the pionic method, the overall systematic error in determining Xqo comes from the un- 
certainty in the IRF including that due to the energy calibration uncertainty (±10%), unknown 
contributions of the overlapping sources (+0/ — 2, +0/ — 19, and +0/ — 8%), uncertainty in the CR 
intensity (±10%), uncertainty in the pp pion production cross-section (±5%), and uncertainty in 
the contribution from heavier nuclei (±5%). We conservatively quote the linear sum of these com- 
binations as the possible systematic error for the 3 Orion regions, which are +30/ — 32, +30/ — 49, 
and +30/ — 38%, as given in Column 5 of Table [2j 

The systematic errors that can affect Xqo differently in the 3 Orion regions are variation in 
the CR intensity within ~ 1 kpc (±5%) and the overlapping sources. The overall error of this kind 
is conservatively estimated to be the linear sum of the two, +5/ — 7, +5/ — 24, and +5/ — 13%. 



4.2. Gamma-ray intensity and E{ J — H) 



The line-of-sight visual attenuation, Ay, are often used as a gas-mass trac er in theory-based 



Burgh et al 



2010 



studies of the CO fraction in all molec ules including carbon and hydrogen (e.g 
Wolfire et al.l |2010| ; I Glover et al.l |2010| . and references therein) . To calibrate crudely our mass col- 
umn density with Ay used in these theory-based analyses, we have related the ga mma-ray count s 
on the horizontal axes of Fig. [4] and E(J — H) in the 3 Orion regions measured by iDobashil (|201ll ). 



We note that the atomic and molecular components are assumed to be contained within a fixed 
length (e.g. 20pc) along the line-of-sight in the theory-based analyses while the components are 
measured as column densities integrated over unknown lengths along the line-of-sight in observa- 
tions. Moreover E{J — H) is known to trace the H2 gas but also pick up some H I gas through dust 
mixed with it. Hence the cross-calibration works at best crudely and only in the regions of clouds 
where the H2 longitudinal distribution is well confined and the H2 volume density dominates over 
that of H I. Despite these uncertainties, it is important that our measurements be compared with 
theory-based analyses. 



We found good linear relations for the pixels with high gamma-ray counts (> 300 per deg 2 ) 
in all 3 Orion regions and could correlate the gamma-ray count scale on the horizontal axes of 
Fig. 0] to E(J — H) assuming Ay = Rv~ejh x E(J — H). T he Rv- e.jh has been deter mined 



observationally and its value ranges between 7.8 (|Dobashill201ll ) to 10.9 (jCardelli et al.lll989l ). The 
highest point in our count map is ~ 700 per pixel in Orion A Region II where H2 concentration 
is highest and the corresponding value of Ay is ~ 5 when averaged over lxl deg 2 pixels for an 
assum ed value of Rv-ejh = 7.8. So Ay = 5 on the horizontal axes of Figs. 5 and 6 in 1 Glover et al.l 
(|2010l ) corresponds crudely to ~ 700 counts per pixel assuming H2 is well confined (e.g., to ~ 20 pc) 
along the line-of-sight. 



4.3. Summary of the Results 

The results obtained in the present work are significant beyond the estimated systematic errors. 
They are: 



1. 



3. 



4. 



Linearity holds between mass density associated with the Orion clouds and Wco' As discussed 
in Subsection [3T2] and shown in Fig. 01 our results suggest th at CRs penetrate t o all translucent 
part of the clouds. Possible shielding of CRs discussed in lAharonianl (|200lh does not apply 
to most parts of the Orion clouds. 

The Xqo factors calculated with the pionic method and with the H2/H I method differ by 
~ 15% but agree within the estimated systematic error (Table [2]). The difference can be 
explained by uncertainties in the column densities of H I and calculation of gamma-ray emis- 
sivity per H I atom. 



The Xqo factor obtained with the H2-template-2 is found to be larger by ~ 40 — 60% in 
Orion A Region I than Orion A Region II and Orion B for the two methods. The difference 
is much larger than the systematic error that can affect the Xco factor differently in the 3 
Orion regions (Table [2]). 

In the "dark gas" scenario, the added "dark gas" accounts for the majority of the gas not 
traced by Wco- One Xqo factor can then describe the Wco-traced H2 distribution in the 
ROI. 



4.4. Interpretation of our results on Xco 

Historically the relation between iV(H2) to Wco has been considered to depend on the envi- 
ronment around the molecular cloud. The environmental factors discussed in the literature are: 



Metallicit y: This possibility h as been discussed in the literature since the late 1980's (e.g.. lElmegreen 
19891 : iBolatto et al.l [1999J) . According to an empirical formula proposed to relate X^q to 



[0/H] (|Wilsonlll995l ; lArimoto et al.lll99q ). the metallicity must be ~ 2 times higher in Orion 
A Region I to account for the observed difference in Xc o between Region I and II, which is 



unlikely according to Galactic-scale measurements (e.g. lEsteban et al.ll2005l ). We note that 
metallicity is generally considered to be an important environmental factor influencing the 
H2-to-H I ratio. 

Overlapping HI clumps: Compact HI clo uds with angular diameters of 1 — 2deg have been 



found in various Galactic locations (e.g. iBraun &: Stromlll986l ; iKavars et al.ll2003l ; iLee et al 



20081 ). A new reanalysis of the LAB H I survey shows no such concentration detected at the 



sensi tivity level of the present study overlapping with the Orion A and B (IKalberla et al 



20101 ). 



Low density H2 not traced well by Wco : Existence of diffuse H2 gas not traced well by Wqo 
has been discussed in the literatu re cited in Section [J an d Fermi analyses are bringing the 
discussion to a quantitative level (jAckermann et al.ll2010l ). We refer to the following recent 
works on the H2 and CO fractions and try to interpret our results: 



Burgh et al.l (|2010l ) have studied the fractions based on Hubble Space Telescope observa- 



tions and characterized the Xco dependency on N(H.2). 



Wolfire et al.l (|2010l ) have studied chemical composition of a model cloud theoretically 



and found that CO becomes depleted because of photodissociation in the periphery 
where the gas density decreases. 



Glover et al.l (|2010l ) have studied the time-dependent H2 and CO fractions in clouds 



through computer simulations and found Xqo increases sharply where N(H.2) decreases 
for Ay < 3.5. 

All of the above studies predict that the CO/(total C) fraction drops as the H2 column density 
decreases, as toward the periphery of Orion A and B. However the Wco-to^ relation and 



the ab undance of H2-without-CO gas may be more complicated. For example, llkeda et al 



(|2002j) found that N(C l)/A r (CO) increases to high values along all of the peripheries whereas 
we find Region I of Orion A to be more abundant in CO-d epleted gas than R egion II. The 
prediction that Xqq increases sharply in regions Ay < 3.5 bv lGlover et al.l (|2010l ) is consistent 
with our finding that the "dark gas" is concentrated in the high-longitude end of Orion A 
where Wqo becomes low. 



5. Conclusion 

We have reported on the first 21 months' observations of Orion A and B with the Fermi 
Gamma-ray Space Telescope in the energy band between ~ 178 MeV and ~ 100 GeV. We have mea- 
sured the mass column density distribution within the clouds at the angular scale of the instrument 
PSF using the pp — > 7 production cross-section accurately calibrated at accelerators as well as using 



the gamma-ray emissivity of the local H I gas. We found with the pionic method that a linear rela- 
tion holds between mass density and W C o with X co = 2.34, 1.43, 1.35 xlO 20 cm" 2 (K km s" 1 )" 1 
with a systematic uncertainty of +5/ — 7, +5/ — 24, and +5/ — 13% (relative in the 3 regions), and 
+30/ - 32, +30/ - 49, and +30/ - 38% (absolute) for Orion A Region I, Region II, and Orion B, 
respectively. These values are consistent with the Xqo values determined with the more traditional 
H 2 /HI method (Xco = 1-97, 1.20, 1.14 xlO 20 cm~ 2 (K km s -1 ) -1 ) within our overall systematic 
error. This implies that Galac tic CRs are penetra ting into most parts of the clouds. The analyses 
also included the "dark gas" (jGrenier et al.l 120051 ) not traced by CO or H I. We found that the 
gamma-ray flux associated with the "dark gas" spatial template exceeds that associated with the 
Wqo template in Orion A Region I. The situation is reversed in Region II and in Orion B. This is 
generally consistent with the fit finding a higher Xco value for Orion Region I in the absence of 
the dark-gas template. 

We have interpreted the increase in Xco an d "dark gas" fraction in Ori on A Region I in th e 
light of recen t stud i es of the relation b etween the H2 and CO fractions by iBurgh et al.l (|2010i ); 
Wolfire et al.1 (|2010l ) ; I Glover et al.l (120101) . Xco is expected to increase rapidly as the gas column 
density decreases to Ay ~ 3.5 or less (jGlover et al.l |2010| ) . The mass column density we have 



measured in Region I corresponds to Ay < 4, close to the predicted threshold for onset of the non- 
linearity predicted between Wco and A^B^). The mass column density drops further (Ay < 2) 
toward the high Galactic longitude end of the Orion A where the gas becomes "dark" to Wqo, 
consistent with the predicted non-linear relation. 

The Fermi-LAT collaboration is continuing to reduce uncertainty in the IRF, identify extended 
gamma-ray sources, and improve the modeling of the Galactic-scale diffuse gamma-ray emission. 
We expect the systematic uncertainties quoted in subsection 14. II to be reduced significantly through 
these efforts. The systematic uncertainty in the CR spectra and the B I mass density also will be 
reduced when the data from new experiments and surveys become available. The present analyses 
can then be updated to a higher precision and the relation among Wco and the gas mass density 
characterized further for various molecular clouds in the Galaxy. 
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Table 1. Gamma-ray emissivity fitted with H2-template-2 



Energy Range (MeV) 



178 


- 237 


237 


- 316 


316 


- 422 


422 


- 562 


562 


- 750 


750 


- 1000 


1000 


- 1334 


1334 


- 1778 


1778 


- 2371 


2371 


- 3162 


3162 


- 4217 


4217 


- 5623 


5623 


- 10000 


10000 


- 23714 


23714 


- 100000 



Emissivity per H a 



4.81 ± 0.26) 


X 


10" 


-29 


3.15 ± 0.10) 


X 


10" 


-29 


1.81 ± 0.04) 


X 


io- 


-29 


1.05 ± 0.02) 


X 


10" 


-29 


5.72 ± 0.12) 


X 


io- 


-30 


3.20 ±0.08) 


X 


io- 


-30 


1.69 ±0.09) 


X 


io- 


-30 


8.75 ± 0.30) 


X 


io- 


-31 


4.19 ±0.25) 


X 


io- 


-31 


1.83 ±0.14) 


X 


io- 


-31 


7.97 ±2.72) 


X 


io- 


-32 


4.07 ±0.29) 


X 


io- 


-32 


1.19 ±0.38) 


X 


io- 


-32 


1.42 ± 1.01) 


X 


io- 


-33 


4.16 ±3.06) 


X 


io- 


-35 



Orion Region I 



1.04 ±0.08) 


X 


io- 


-28 


6.36 ±0.39) 


X 


10" 


-29 


3.68 ±0.21) 


X 


io- 


-29 


1.95 ±0.11) 


X 


io- 


-29 


1.29 ±0.07) 


X 


io- 


-29 


5.97 ±0.37) 


X 


io- 


-30 


3.16 ±0.23) 


X 


io- 


-30 


1.55 ±0.13) 


X 


io- 


-30 


7.49 ±0.77) 


X 


io- 


-31 


4.01 ±0.47) 


X 


io- 


-31 


2.16 ±0.29) 


X 


io- 


-31 


6.98 ± 1.46) 


X 


io- 


-32 


2.27 ±0.47) 


X 


io- 


-32 


2.12 ±0.95) 


X 


io- 


-33 


1.07 ±0.79) 


X 


io- 


-34 



Emissivity per Wco b 
Orion Region II 



(5.20 ±0.84) 


X 


io- 


29 


(3.50 ±0.37) 


X 


io- 


-29 


(2.16 ±0.18) 


X 


io- 


29 


(1.24 ±0.09) 


X 


io- 


-29 


(6.15 ±0.46) 


X 


io- 


-30 


(4.08 ±0.27) 


X 


io- 


-30 


(2.08 ±0.15) 


X 


io- 


-30 


(1.06 ±0.08) 


X 


io- 


-30 


(6.08 ±0.53) 


X 


io- 


-31 


(2.60 ±0.29) 


X 


io- 


31 


(1.23 ±0.17) 


X 


io- 


31 


(4.77 ±0.95) 


X 


io- 


-32 


(1.05 ±0.27) 


X 


io- 


32 


(1.48 ±0.53) 


X 


io- 


-33 


1.57 ±3.34) 


X 


io- 


-34 



Orion B 



6.17 ± 0.52) 


X 


10 


29 


3.64 ±0.26) 


X 


io- 


-29 


2.32 ±0.13) 


X 


io- 


-29 


1.17 ±0.07) 


X 


io- 


-29 


6.96 ±0.41) 


X 


io- 


-30 


3.50 ±0.24) 


X 


io- 


-30 


1.70 ±0.14) 


X 


io- 


-30 


8.71 ±0.80) 


X 


io- 


-31 


4.82 ± 0.49) 


X 


io- 


-31 


1.92 ±0.27) 


X 


io- 


-31 


9.34 ± 1.56) 


X 


io- 


-32 


4.33 ± 0.90) 


X 


io- 


-32 


1.10 ±0.28) 


X 


io- 


-32 


1.61 ±0.59) 


X 


io- 


-33 


1.14 ±0.53) 


X 


io- 


-34 



Note. — Errors are statistical only 
a MeV~ 1 s~ 1 sr~ 1 per H atom 

b MeV- 1 s- 1 sr- 1 (2 X 10 20 c m - 2 (K km s" 1 )- 1 )- 1 



Table 2 


. Xqq, obtained 


on H2- 


template- 1, 


2, and 3 






Region 


X CO a on B/2A 


Sys. error 6 (%) 


Xco a on pion 


Sys. error c (%) 


H2-template-l 


Entire ROI 


1.36 ± 0.02 stat 


NA 




1.63 ± 0.02 stat 


NA 




H2-template-2 


Orion A Region I 


1.97 ± 0.05 stat 


+25/ 


- 28 


2.34 ± 0.05 stat 


+30/ 


- 32 


Orion A Region II 


1.20 ± 0.03 stat 


+25/ 


-44 


1.43±0.04 stat 


+30/ 


-49 


Orion B 


1.14 ± 0.03 stat 


+25/ 


- 33 


1.35 ± 0.03 stat 


+30/ 


- 38 


Elsewhere 


1.43 ± 0.04 stat 


NA C 




1.69 ± 0.04 stat 


NA d 




H2-template-3 


Entire ROI 


1.21 ± 0.02 stat 


+25/ 


- 37% e 


1.32 ± 0.02 stat 


+30/ 


-40 e 



a In unit of 10 20 cm~ 2 (K km s r ) x . 

b The systematic error is discussed in Subsection 14.11 it comes from a combination of 
uncertainties in the H I spin temperature and in the fitting process. The systematic errors 
which may apply differently to the 3 Orion regions are +5/ — 8, +5/ — 24, and +5/ — 13%, 
respectively. 

c The systematic error is discussed in Subsection 14.11 The systematic errors are the same 
as b. 

d We have not attempted to estimate systematic error outside of the Orion regions in 
this study. 

e The average of the systematic errors estimated for the 3 Orion regions. 



Table 3. Gamma-ray emissivity fitted with H2-template-3 



Energy Range (MeV) Emissivity per H atom a Emissivity per Wco b Emissivity per E(B — V)res c 



178 


- 237 


(4.51 ±0.08) 


X 


10" 


-29 


(5.56 ±0.31) 


X 


10" 


-29 


(1.00 ± 0.10) 


X 


10" 


-27 


237 


- 316 


(2.99 ±0.08) 


X 


io- 


29 


(3.39 ±0.27) 


X 


10" 


-29 


(5.60 ± 0.83) 


X 


10" 


-28 


316 


- 422 


(1.68 ±0.07) 


X 


io- 


29 


(2.03 ±0.11) 


X 


10" 


-29 


(3.95 ± 0.36) 


X 


io- 


-28 


422 


- 562 


(1.02 ±0.08) 


X 


10" 


29 


(1.13 ±0.05) 


X 


10" 


-29 


(2.07 ±0.16) 


X 


io- 


-28 


562 


- 750 


(5.39 ±0.08) 


X 


10" 


30 


(6.51 ±0.21) 


X 


10" 


-30 


(1.37 ±0.09) 


X 


io- 


-28 


750 - 


- 1000 


(2.97 ±0.09) 


X 


10" 


-30 


(3.54 ±0.16) 


X 


10" 


-30 


(6.57 ±0.60) 


X 


io- 


-29 


1000 


- 1334 


(1.58 ±0.05) 


X 


10" 


-30 


(1.86 ±0.09) 


X 


10" 


-30 


(3.57 ±0.35) 


X 


io- 


-29 


1334 


- 1778 


(8.00 ± 1.02) 


X 


10" 


-31 


(9.37 ± 0.43) 


X 


10" 


-31 


(1.86 ± 0.16) 


X 


io- 


-29 


1778 


- 2371 


(3.64 ±0.25) 


X 


10" 


-31 


(5.00 ±0.31) 


X 


10" 


-31 


(7.45 ± 1.19) 


X 


10" 


-30 


2371 


- 3162 


(1.51 ±0.14) 


X 


10" 


-31 


(2.19 ±0.17) 


X 


10" 


-31 


(4.82 ± 0.68) 


X 


io- 


-30 


3162 


- 4217 


(6.56 ±0.89) 


X 


10" 


-32 


(1.06 ±0.10) 


X 


10" 


-31 


(2.18 ± 0.40) 


X 


io- 


-30 


4217 


- 5623 


(3.82 ± 1.71) 


X 


10" 


-32 


(4.31 ± 0.49) 


X 


10" 


-32 


(6.50 ± 2.34) 


X 


10" 


-31 


5623 - 


- 10000 


(1.06 ±0.14) 


X 


10" 


-32 


(1.07 ±0.16) 


X 


10" 


-32 


(2.20 ± 0.68) 


X 


10" 


-31 


10000 


- 23714 


(1.35 ±0.15) 


X 


10" 


-33 


(1.68 ±0.24) 


X 


10" 


-33 


(1.72 ± 0.92) 


X 


10" 


-32 


23714 


- 100000 


(4.62 ±6.52) 


X 


10" 


-35 


(9.55 ±3.50) 


X 


10" 


-35 


(1.46 ± 1.22) 


X 


io- 


-33 



Note. — Errors are statistical only 
a MeV -1 s~ 1 sr~ 1 per H atom 

b MeV"" 1 s- 1 sr- 1 (2 X 10 20 cm- 2 (K km s" 1 )- 1 )- 1 
c MeV~ 1 s- 1 sr- 1 (2 X lO^mag) -1 



Fig. 1. — Gamma-ray count distribution in the Orion region in the energy band between 178 MeV 
and 100 GeV in the Hammer-Aitoff projection on the Galactic coordinates. The pixel size is 
0.2 x 0.2 deg 2 . 
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Fig. 2. — (a) iV(H I) template summed over the line-of-sight velocity. The pixel size is 0.5° x 0.5°. 
The dashed lines show the boundary of the ROI. (b) Wq q template used in H2-template-l and 
H2-template- 2. We used NANTE N data (jFukui et al.ll201ll ) in the area bounded by the solid lines 
and those by bame et all (|200ll ) elsewhere. Pixel resolution is 0.125° x 0.125°. (c) E{B - V) Tes 
template used in H2-template-3. Pixel resolution is 0.5° x 0.5°. 
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Fig. 3. — (a) Gamma-ray count distribution in the ROI after subtracting the fitted "background" 
that is the sum of the H I, IC, point-source, and isotropic components, (b) The fitted model 
map obtained by assuming one common Xco for the ROI (H2-template-l). Dashed lines define 
the boundaries of the 3 Orion regions, Orion A Region I, Region II and Orion B. (c) Same as 
(b) but obtained by assuming 4 different Xqo f° r Orion A Region I, Region II, Orion B, and 
elsewhere (H2-template-2). (d) Same as (b) but obtained by adding E(B — V) Tes to BVtemplate-l 
(H2-template-3). 




Fig. 4. — (a) Correlation between the gamma-ray count distribution shown in Fig. [JJt (the hor- 
izontal axis) and that fitted with H2-template-l in Fig. [3b (the vertical axis) for all pixels in the 
3 Orion regions. Points represent pixels in Orion A Region I (black circles), Region II (red tri- 
angle), and Orion B (blue squares) with fitted lines black, red, and blue, respectively. Error bars 
represent statistical errors in counts in pixels. Same after replacing the vertical axis for that fitted 
with H2-template-2 (b) and for that fitted with H2-template-3 (the sum of Wqo an d E(B — V) Tes 
components) (c). 
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Fig. 5. — Difference in log(Likelihood) between the spatial fit using H2-template-3 (dotted line) 
and either that with H2-template-l (dashed line) or that with H2-template-2 (solid line) in the ROI 
for the 22 energy bins. Note that the lines are drawn between the data points only to guide the 
eye. 




Fig. 6. — Residue in the energy-summed gamma-ray counts of the spatial fit with H2-template-l 
(a), H2-template-2 (b), and H2-template-3 (c), binned in 1 x 1 deg 2 pixels. The black dotted lines 
show the boundaries of the 3 regions, Orion A Region I, II, and Orion B. 




Fig. 7. — Spectral energy densities (SED) associated with local H I (Ts = 125 K assumed) (a) and 
that associated with H2-template-l (b). The lines are: total (solid), bremsstrahlung (dotted) and 
pion decay (dashed). The CR spectral shape and electron-to-proton ratio at the Orion clouds were 
fixed to those used by GALPROP. The vertical axes are normalized to the column density of H I in 
unit of 1cm -1 for (a) and to 2 x Xco m unit of 10 20 cm _2 (K km s -1 ) -1 for (b). The energy bins 
between No. 13 and No. 22 are combined to wider energy bins. Vertical bars represent statistical 
errors. Note that the spectral fit to H I is not used in evaluating Xco- 
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Fig. 8. — Spectral energy density (SED) associated with local H I (T$ = 125 K assumed) (a), Orion 
A Region I (b), Region II (c), and Orion B (d) fo r the fit with H2-template-2. The corresponding 
SED obtained for the local H I (jAbdo et al.ll2009d ) is shown by gray squares in (a). The assumption 
about the CR, the line legends, and the vertical axis units are the same as in Fig. [7J 




Energy (MeV) 



Fig. 9. — SED associated with local H I (Tg = 125 K assumed) (a), that associated with Wco 
(b), and that associated with E{B — V) res (c) obtained with H2-template-3. The line legends and 
vertical axis units are the same as in Fig. [7J 
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Fig. 10. — Gamma-ray spectra spatially associated with two H2 templates in the 3 Orion regions 
marked in Fig. [3b: (a) the sum of the 3 regions obtained with H2-template-2; (b) the sum of the 3 
regions with H2-template-3; (c) Orion A Region I obtained with H2-template-3. Black circles show 
the isotropic component, red squares H I, green upward triangles CO, and purple dashed line the 
inverse Compton. Blue downward triangles in (b) and (c) represent the spectra associated with 
E(B - V) ms . 



